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A WEIGHTED MINIMUM GRADIENT PROBLEM WITH 
COMPLETE ELECTRODE MODEL BOUNDARY CONDITIONS FOR 
CONDUCTIVITY IMAGING* * * § 

ADRIAN NACHMANt, ALEXANDRU TAMASAN*, AND JOHANN VERAS§ 


Abstract. We consider the inverse problem of recovering an isotropic electrical conductivity 
from interior knowledge of the magnitude of one current density field generated by applying current 
on a set of electrodes. The required interior data can be obtained by means of MRI measurements. On 
the boundary we only require knowledge of the electrodes, their impedances, and the corresponding 
average input currents. From the mathematical point of view, this practical question leads us to 
consider a new weighted minimum gradient problem for functions satisfying the boundary conditions 
coming from the Complete Electrode Model (CEM) of Somersalo, Cheney and Isaacson. We show 
that this variational problem has non-unique solutions. The surprising discovery is that the physical 
data is still sufficient to determine the geometry of (the connected components of) the level sets 
of the minimizers. We thus obtain an interesting phase retrieval result: knowledge of the input 
current at the boundary allows determination of the full current vector field from its magnitude. We 
characterize locally the non-uniqueness in the variational problem. In two and three dimensions we 
also show that additional measurements of the voltage potential along a curve joining the electrodes 
yield unique determination of the conductivity. The proofs involve a maximum principle and a new 
regularity up to the boundary result for the CEM boundary conditions. A nonlinear algorithm is 
proposed and implemented to illustrate the theoretical results. 


Key words, minimum gradient, conductivity imaging, complete electrode model, current den¬ 
sity impedance imaging, minimal surfaces, magnetic resonance electrical impedance tomography, 
current density impedance imaging 
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1. Introduction. We consider the inverse problem of reconstructing an inho- 
mogenuous isotropic electrical conductivity a in a domain C K™, n > 2, from 
interior knowledge of the magnitude a of one current density field and of correspond¬ 
ing boundary data. 

Most of the existing results on this problem (see a brief survey of previous work at 
the end of this introduction) consider Dirichlet boundary conditions. In this paper we 
study boundary conditions which model what can actually be measured in practical 
experiments. We work with the beautiful Complete Electrode Model (CEM) originally 
introduced in [25] and shown to best describe the physical data: For k = 0,..., N, let 
e/c C dtt denote the surface electrode of constant impedance Zk through which one 
injects a net current Ik- 

The CEM assumes the voltage potential u inside and the constant voltages Uk s 
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on the surface of the electrodes distribute according to the boundary value problem 
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where v is the outer unit normal. If a solution exists, an integration of EH over 
together with EH and EH show that 

N 

E J * = o (i-5) 

k—0 

is necessary. Physically, the zero sum of the boundary currents account for the absence 
of sources/sinks of charges. The constants Uk appearing in (11.211 represent unknown 
voltages on the surface of the electrodes, and the difference from the traces u\ ek of 
the interior voltage potential governs the flux of the current through the skin to the 
electrode. We refer to the problem EH, EH, EH, and EH as the forward problem. 

Under the assumptions that Q is a Lipschitz domain, the conductivity is essentially 
bounded with real part bounded away from zero, the electrodes e*, are (relatively) open 
connected subsets of dfl whose closure are disjoint, the impedances Zk have positive 
real part, and the injected currents Ik satisfy EH, the forward problem has a unique 
solution (u; (Uq, ■■■,Un)) £ 7? 1 (U) x C^ -1 " 1 up to an additive constant, as shown in 
[28] . We normalize this constant by imposing the electrode voltages U = {Uq, ..., Un) 
to lie in the hyperplane 

n := {U £ K jV+1 : U 0 + ... + U N = 0}. (1.6) 

The net input currents Ik, k = 0, ...,7V as in EH, generate a current density field 
J = — erVu, where (u, U) £ 7V 1 (U) x II is the solution of the forward problem. 

In this paper we consider the inverse problem of determining er, given the magni¬ 
tude 


a = |crVrt| 


(1.7) 


of the current density field inside Q. 

The conductivity a is unknown but assumed real valued and satisfying 

essinfQcr(:r) > 0. (1.8) 

Each electrode ek C dfl, k = 0, ...,7V, is a known Lipschitz domain subset of the 
boundary. The surface impedances Zk are assumed real valued. In general we allow 
them to be inhomogenous functions on the electrodes satisfying 

essinf efc 2fc > 0, for A; = 0,..., TV. (1.9) 

Further smoothness conditions will be assumed for some of the results in this paper. 
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We note that, in practice, interior measurements of all three components of the 
current density J can be obtained from three magnetic resonance scans involving two 
rotations of the object m- However recent engineering advances in ultra-low field 
magnetic resonance may be used to recover J without rotating the object [26]. We 
hope that the results presented here may lead to further experimental progress on 
easier ways to measure directly just the magnitude of the current. 

We start by remarking that there is non-uniqueness in the inverse problem stated 
above, as can be seen in the following example: Let H = (0,1) x (0,1) be the unit 
square. We inject the current I\ = 1 through the top electrode ei = {(x, 1) : 0 < x < 
1} of impedance Z\ > 0, “extract” the current Iq = —1 through the bottom electrode 
e o = {( 0 ,x) : 0 < x < 1} of impedance Zq = Z\ + 1, and measure the magnitude 
a = 1 of the current density field in S2. Then, for every ip : [0,1] —> [y>(0), y>(l)] 
an increasing Lipschitz continuous function, satisfying <p(0) + y>(l) = 1, the function 
u v (x,y) := ip(y) solves the forward problem (11.11) . (11.21) . (11.31) . and dm corresponding 
to a conductivity a v (x, y) = 1 /<p'(y), yet the magnitudes of the corresponding current 
densities yield the same interior measurements ojVu| = cr^Vu^l = 1. 

More generally, if (u, U ) £ H 1 ^) x n is the solution of the forward problem for 
some cr, let ip £ Lip(u(Cl )) be any Lipscliitz-continuous increasing function of one 
variable, such that <p(t) = t + Ck whenever t £ u(ek), for each k = 0 and 
constants Ck satisfying Ylk =o Ck = 0. One can easily verify that the function 


( 1 . 10 ) 


u v := ip o u 


solves the forward problem with the conductivity 


while ojVxt| = ij^IVu^I. 

For Holder-continuous conductivities, in Theorenr l3.1l we prove that (11.101) . (11.111) 
must hold in a neighborhood of any non-critical point. The following example shows 
that ( 11 . 101 ) . ( 11 . 111 ) need not hold in the whole domain Q for a single function <p. 

Let SI C U 2 be the curvilinear octagon with hyperbolic sides obtained from the 
unit disc, by carving out the peripheral regions along the branches of the hyperbolas 
x 2 — y 2 = ± 7 j, and xy = ±-^. The electrodes eo,e 2 are defined respectively by 
each connected component of x 2 — y 2 = | inside the disc, whereas ei,e 3 are defined 
respectively by the connected components of x 2 — y 2 = Thus defined, all the 
electrodes have equal length, denoted by l. We assume constant impedances Zq = 
zi = Z 2 = Z 3 = 1. Through eo and e 2 one inputs the net currents Iq = I 2 = an d 
extracts Ii = I3 = —d through ei and e 3 . In S2 one measures \J(x,y)\ = 2^/x 2 + y 2 . 
It is easy to check that the constant conductivity a = 1 is a possible solution to the 
inverse problem, with the corresponding voltage u(x,y) = x 2 — y 2 . However, this is 
not the only possibility: Let 


w± := {(x,y) £ 12 : e < x 2 - y 2 < i - e, ± x > 0 }, 


and 1 p± : (0, —> (0, be any two increasing Lipschitz continuous functions with 
<P±{t) = t in (0, e) U (i — e, for some e > 0 sufficiently small. Define a new 

conductivity by 



( 1 . 12 ) 
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As in the previous examples, one can check that a is also a solution of the inverse 
problem for the same \J(x,y)\ = 2^/x 2 + y 2 1 with corresponding potential v equal to 
u on Q \ ui±, and equal to (p± o u on ui±, respectively. 

Similar to the approach in [22] . we formulate the inverse problem in terms of a 
weighted minimum gradient problem. Here we need a functional which is appropriate 
for the boundary conditions coming from the Complete Electrode Model and is entirely 
defined in terms of the data in the inverse problem. To obtain such a functional, 
we found it necessary to revisit the forward problem and recast it as a minimization 
problem, see Appendix A. We are then able to show that the solution (it, U) £ H 1 (f2) x 
n of the forward problem is a global minimizer of the functional 


N N 

G a (v,V)= / a\Vv\dx + Y. / ^(v-V k fds-Y,hV k 

fc=0-' efc " k fc=0 


(1.13) 


over H 1 (H) x n, with a = <t|Vu| as in (11.71) . 

Using G a , we found the surprising fact that, given the positions and impedances 
of the electrodes, knowledge of the magnitude of one current density field and of 
the corresponding average applied currents (just one number in the case of two elec¬ 
trodes!) is still sufficient to determine the geometry of the connected components of 
the equipotential sets. Furthermore, we remark that, since we recover the direction 
(including orientation) of the electric field Vit, we also obtain an interesting phase 
retrieval result: the full current density vector field J = is recovered from its 
magnitude a, and knowledge of the input currents Ik on the surface electrodes, even 
though the conductivity is not uniquely determined. 

Uniqueness of the conductivity can be restored by additional measurement of 
the voltage potential on a curve T connecting the electrodes; see Theorem 14.31 The 
additional measurement involves only a one dimensional subset of boundary measure¬ 
ments, much less than required in the existing results for the Dirichlet problem. 

For the unique determination result we assume that the curve T and the set of 
electrodes satisfy a topological assumption, see EH)- 

To illustrate the theoretical results, we propose an iterative algorithm and perform 
a numerical experiment, see Section [5] Similar to the algorithm in |22| we decrease 
the functional G a on a sequence of solutions of forward problems for updated conduc¬ 
tivities. 

Conductivity imaging using the interior knowledge of the magnitude of cur¬ 
rent densities was first introduced in m- The examples of non-existence and non¬ 
uniqueness for the ensuing Neumann problem lead the authors of [12] to consider the 
magnitudes of two currents. The possibility of conductivity imaging via the mag¬ 
nitude of just one current density field was shown in |21j via the Cauchy problem, 
and in [22] 124] via a minimum gradient problem with Dirichlet boundary conditions. 
Existence and uniqueness of such weighted gradient problems was studied in [8] . Ex¬ 
tensions to the case of inclusions with zero or infinite conductivity were obtained in 
ESI EDI- a structural stability result for the minimization problem can be found in 
[25] . Reconstruction algorithms based on the minimization problem were proposed 
in [22] and m, and based on level set methods in EH [2! [30]. A local Holder- con¬ 
tinuous dependence of er on |J| (for unperturbed Dirichlet data) has been recently 
established in m- For further references on determining the isotropic conductivity 
based on measurements of current densities see eh m m m hh m uni , and for 
some results on anisotropic conductivities see MEH- 



A weighted minimum gradient problem in conductivity imaging 


5 


2. A weighted minimum gradient problem for the CEM boundary con¬ 
ditions. In this section we show that the solution of the forward problem is a global 
minimizer of the functional G a in (11.131) over U 1 (fl) x II. The regularity assumptions 
are the ones from the forward problem. 

Proposition 2.1. Let f1 C n > 2, be a bounded domain with Lipschitz 
boundary. Let ek, for k = 0, ...,1V, be disjoint subsets of the boundary of positive 
surface measure, a £ L°°(f2) and Zk £ L°°(ek) be bounded away from zero, and Ik 
satisfy (EH)- Let ( u , U) £ H l (Lt) x II be the unique solution of the forward problem 
(11.11) . (11.21) . (11.31) . (11.41) . If a := cr|Vu|, then 

G a (v, V) > G a (u, U), V(v, V) £ H\n) x II. (2.1) 


Proof. For any ( v , V ) £ H 1 ( f2) x II, we have the inequality 


1 N N 

G a {v,V)= / a\\7u\\Vv\dx + - / — {v - V k ) 2 ds - Y' I k V k 

Jn t'o 
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> / crVzt • Vvdx + / 

fc _ Q Je k 


_L 


ds 


f du u V' 

./an ^ 


T ( » - v t f - 
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= E 


fc =0 efc 
AT 


fc =0 efe L 
<9u 1 . Tr , 2 Tr du 

m TPirj v - Vt) - v “ a di 


ds 


ds 


= £ 


fc =0 l/efc 


1 

Zfc 


1 


(w - Vk)(Uk -u) + ^(v - Vfc) 2 


ds, 


( 2 . 2 ) 


where the first equality uses s the next line uses ESI and the Cauchy-Schwarz 
inequality, the next equality uses ed and the divergence theorem, the third equality 
uses El, and the last equality uses E^- 

In particular, when v = u, the inequality in the estimate (12.21) holds with equality 
yielding 


and 


N 


G a (u,V) = J2 


k=0 ek 


Zk 


-(U - Vk ) 2 - {u - V k )(u - Uk) 


ds, 


(2.3) 


Ga(u, U) 



Ukfds. 


(2.4) 


The global minimizing property m then follows from (12.21) and (12.41) using the 
pointwise inequality 

^(v - V k ) 2 - (v- 14)(u - U k ) >-^(u- U k ) 2 . (2.5) 

□ 

Note how each of the CEM equations EH , El , ESI , EH was used in the 
above proof. 
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3. Local characterization of non-uniqueness and applications. In this 
section we state and prove our main result and its consequences to the conductivity 
imaging problem. 

Theorem 3.1. Let S2 c d > 2 be a bounded, Lipschitz domain, and let e k , 
k = 0, ...,1V, be disjoint subsets of the boundary of positive surface measure. Assume 
that the corresponding impedances z k satisfy (ED, and that the given currents Ik 
are such that EU holds. Let (u,U),(v,V) £ H 1 (fl) x II, be the solutions of the 
forward problem EEJ. E2D- m, and E3D corresponding to unknown conductivities 
a, a £ C a (Ll) satisfying (11.81) . Assume that 

<t|Vu| = ct|Vu| > 0 a.e.in S7. (3-1) 


Then: 

(i) for each k = 0 

u\e k - u k = v\e k - Vk, a.e. on e k - (3.2) 

(ii) for a.e. x° £ S2 with |Vu(a;o)| ^ 0, there exists a neighborhood O o of x° and 
a function ip £ C 1 (u(Otf)), such that 

v = ip o u, in Oq, (3-3) 

and 

a = a , in O 0 . (3.4) 

p' o u 

Moreover, the equalities (ED and ED hold on the union of the connected components 
of the level sets of u passing through Oq. 

Proof. From the interior elliptic regularity we have u,v £ C 1,a (Ll); see, e.g., 
[4] Theorem 8.34]. Thus, the sets of critical points of u and v are closed, and, by 
hypothesis (13.11) they are negligible. Let S denote their union. It follows that H\ S 
is open and dense in SI. 

According to Proposition 12.11 (u.U) and (v,V) are both minimizers of G a , and 

thus 


G a (u,U) = G a ( v,V). (3.5) 

In particular, using (12.41) . (13.51) . (12.21) . and (12.51) . the inequalities 


1 


N 


1 


—(u — U k ) 2 ds = G a (v, V) 
' 1 Zk 


k=0 Jek 


N r i 

i„ _r» J eu * k 


k=0 Jek 
N 


('V - V k )(U k — u) + -{v - 14 ) 2 


ds 


> -9 / —( u-Ukfds 


fc= o Jek 


Zk 


must be equalities. Thus, 

N 


Vk ) - ( w - u k )} 2 ds = 0. 

7_n J eu ZZk 


k—0 Jefc 


(3.6) 





A weighted minimum gradient problem in conductivity imaging 


7 


Now (13.61) together with Zk > 0 yield (13.21) for each k = 0,1, N. 

We will prove part (ii) for any point in f l \ S. 

Since m holds with equality, we must also have 

Vti • = |Vm| • |Vv|, a.e. on f2. (3.7) 

Since both gradients are continuous, in view of m they must be parallel whenever 
one of them is nonzero, in particular 


Vu = ;iVu in \ S, 


(3.8) 


for some y. £ C(fl \ S) with /z > 0 in Cl \ S. Moreover, since 0 \ S' is dense, y extends 
by continuity to the whole domain Cl. 

Let x° £ Cl \ S be arbitrarily fixed. For some component Xj, t^-(x°) > 0. 

Consider G = (G 1 , ...,G n ) : 12 —> R" defined by G k (xi ,..., x n ) = Xk if k ^ j, and 
Gi{x i,... ,x n ) = u[x i, ...,x n ). Then the Jacobian determinant 


\DG(x°)\ = 




> 0 , 


and, by the inverse function theorem, there is a neighborhood Oq of a"o in \ S, such 
that G : Oq — > G{Oq) is a diffeomorphisnr. 

Consider v o G _1 : G(O 0 ) —> K. We use (13.81) to show that 

G(Oo) 9 y = (yi,.... ,y n ) 9so G~ 1 (y) 

is independent of yk, for all k ^ j. Indeed, for y £ G{Oq), 


= o G^ 1 )( y ) 

= /x(G^ a (y))|^ = 0, for k ± j, 

where the second equality uses (13.81) and the next to the last equality uses the definition 
of G. 

For each yj £ v{Oq) we can now well define 


:=voG 1 ( 2 / 1 , y n )- 

If x = G~ 1 (y) £ O 0 , then y.j = G^(x) = u(x), and the equation above shows 


(3.9) 


ip(u(x)) = v(x), for all x £ Oq. 


By @3D we also obtain 

<r(:r)| Vu(ir)| = cr{x)\ Vu(a;)| = a(x)ip' (u(x))\ Vit(x)|, \/x £ Oq. 
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Since |Vu| > 0 in O 0 , the relation (13.41) follows. 

Finally, since u and v are constant on each other’s connected components of level 
sets within f l\S (as can be seen from a differentiation in the direction tangential to the 
level set and (13.81) 1. the identity (13.31) extends to points on any connected component 
of a level set of u passing through Oq. 

□ 

The result above implies that knowledge of the input currents at the boundary is 
sufficient to determine the full current density from measurements of its magnitude 
in the interior, even when the conductivity is not determined uniquely. Thus we have 
obtained the following “phase retrieval” result, which may be of independent interest. 

Corollary 3.2 (Phase retrieval). Let 14 c d > 2 be a bounded, Lipschitz 
domain, and let ek, k = 0, ...,1V, be disjoint subsets of the boundary of positive sur¬ 
face measure. Assume that the corresponding impedances zk satisfy (11.91) . and that 
the given currents Ik are such that (11.51) holds. Let (u,U), (v,V) £ H 1 (Ll) x II be 
the solutions of the forward problem m , & (Ei: and m corresponding to 
unknown conductivities a,d £ C a (Ll) satisfying (11.81) . Let J := aVu and J := &S7v 
be the corresponding current densities. If 

| J| = | J| > 0 a.e.in 14, (3.10) 

then 

J = J in LI. (3.11) 


Proof. For a; 0 £ 14 \ S arbitrarily fixed, let O o and p : v(Oq) —>■ R. be the function 
provided by Theorem 13.11 From (13.41) and (13.31) we have for all x £ Oq, 

J{x) = a(x)S7v(x) =-— Vv(x) = <j(x)\7u(x) = J(x). 

ip' o u{x) 

Since x° is arbitrary in 14 \ S', J and J are continuous in 14, and 14 \ S is dense in 14, 
the result follows. 

□ 

Remark: Theorem !!. H eart (ii), and Corollarv l3.2l are also valid for the Neumann 
problem. The proofs are essentially the same as above, using instead the functional 
G^ ( v ) = f n a\\7v\dx — f gn fvds, where / = is the imposed current on the 

boundary. 

4. Unique determination for two and three dimensional models. To de¬ 
termine the conductivity uniquely, we will assume additional knowledge of the voltage 
potential on some part of the boundary. For brevity, let E := (J^lo ek denote the set 
of all electrodes. 

We show below that knowledge of the voltage potential u along a boundary curve 
r C dLl\E which joins the electrodes is sufficient to yield uniqueness in two and three 
dimensional domains. We assume that T satisfies the topological assumption: 

ruU is a connected set and each connected component of dLl\(TUE) is simply connected. 

(4T) 

In two dimensions the second assumption is trivially satisfied. 

For the proof of the uniqueness result we need continuity of solutions up to the 
boundary for cr-harmonic functions satisfying CEM boundary conditions. 
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As a direct consequence of the Proposition IB. II and the Sobolev embedding the¬ 
orem we have the following. 

Corollary 4.1. Let LI C R” be a C 2 -domain, n = 2,3, and let a £ C a (fl) be 
C 2 -smooth near the boundary. Assume that the electrodes e& has Lipschitz boundaries, 
and Zk £ C 2 (ek), k = 0,..., N. Let (it, U) £ H 1 (LV) x II be the solution to the forward 
problem. Then 

a) u £ C a (Ll), for 0 < a < 1/2. 

b) For any xo £ dLl \ dE, there is a neighborhood Vo C LI of x o, such that 
u £ C 1,01 (Vq), for 0 < a < 1/2. 

Another idea in the uniqueness result below is that the range of it on the union of 
the curve T and the electrodes is the same as the range of u in LI. This follows from 
the following maximum principle for the CEM, which may be of independent interest. 

Proposition 4.2. (Maximum principle for CEM) Under the smoothness as¬ 
sumptions in the Corollary the solution u achieves its minimum and maximum 
on E. I/r C dLt\E is a curve such that TUE is connected, then the range of u over 
ru£ coincides with the range of u over LI. 

Proof. By the weak maximum principle, the maximum M and minimum m of it 
over LI occur on the boundary. By Hopf’s strong maximum principle, at a point of 
maximum, say Xq £ dLl, the normal derivative ^(xo) must be strictly positive. From 
the boundary condition (El we then deduce Xo £ E. The same argument applies to 
a point of minimum, where the normal derivative is strictly negative. Since r U E is 
connected, and u is continuous on LI, it(PU-E) is a closed interval. Since the maximum 
M and minimum m occur on E, then the range u(T U E) = [m, M\ = u(Ll). 

□ 

We are now ready to prove our main uniqueness result. 

Theorem 4.3 (Unique determination). Let U C R" be a C 2 -domain, n = 2,3. 
Assume that the electrodes have Lipschitz boundaries, and Zk £ C 2 (ek), k = 0 ,...,N. 
For the currents Iq,...,In satisfying (11.51) . let (u,U),(v,V) £ iJ 1 (U) x II be the so¬ 
lutions of the forward problem DHL m, HU and m corresponding to unknown 
conductivities a, a £ C a (Ll), which are assumed C 2 -smooth near the boundary and 
satisfying » 

If 

<r|Vtt| = d|Vti| > 0, a.e. in LI, and (4-2) 

u|r = u|r, (4.3) 

then 

u = u in LI, (4-4) 

a = a in LI. (4-5) 


Proof. We will give the proof for the three dimensional case and indicate where 
arguments simplify in the two dimensional case. 

From Corollary 14.11 we know that u and u are continuous up to the boundary. In 
particular, the identity (|3.2I) in Theorem 13.II shows that u = u + Ck on each electrode 
e/c, k = 0,..., N. Since u = u on T D dE by hypothesis (14.31) . we conclude that Ck = 0 
for each k = 0,..., N. So far we showed that it and u coincide on T U E, and following 
Proposition H21 u(Ll) = u(Ll) =: [ m,M ]. 
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We refer to a value t £ [m, M] as being regular if the corresponding t-level set is 
free of singular points. For t-regular, let £ t be a connected component of the t- level 
set. The arguments in [22], Theorem 1.3] showing that £ t reaches dtt do not use any 
boundary information, and thus they remain valid for the CEM boundary conditions; 
we recall them in Appendix B for completeness. 

To prove unique determination, it now suffices to show that Et intersects TUF. 
We reason by contradiction: Assume that £* misses TUE Then the intersection 
Et D dfl is entirely contained in a connected component O C <9f2 \ (HJE 1 ). By 
hypothesis ED O is simply connected. Moreover, as transversal (in fact orthogonal) 
intersection of C' 1 -smooth surfaces, the set E t fldf! is a one dimensional immersed C 1 - 
submanifold without boundary, i.e. a closed curve (in two dimensions it consists of 
two points). Since £ t has no singular points, the curve Etndfl has no self-intersection 
and thus is a simple closed curve embedded in the simply connected subset O. By 
the Jordan curve theorem, £ t D d£l separates O in two parts, one of which, say 0 + , 
is enclosed by £ t D dCl. Let Q + be the subset of 0 whose boundary is E t (J 0 + , and 
define the new function 


-M=={f ) ' ^{LY 4, 

Note that ut may have modified values at the boundary, but only off the electrodes. 
Since f2+ is an extension domain (9f 2+ D Cl = E t has a unit normal everywhere) 
the new map ut £ H 1 (Cl) and strictly decreases the functional (11.1311 unless u = u t . 
This contradicts the minimizing property of u. Therefore u = t in fi + , which now 
contradicts ED- Therefore T, t intersects T(J E, and thus u\s t = u\s t = t. Since the 
set {£ t : t — regular value} is dense in 0, the identity (14.41) follows. Now (14.21) yields 
that a = d, a.e. in f2, and by continuity in f2. 

□ 


5. A minimization algorithm for the weighted gradient functional with 
CEM boundary constraints. In this section we propose an iterative algorithm 
which minimizes the functional G a in (11.1311 . It is the analogue of an algorithm in 
[22l adapted to the CEM boundary conditions. 

The following lemma is key to constructing a minimizing sequence for the func¬ 
tional G a . 

Lemma 5.1. Assume that v £ H 1 ^) satisfies 


e < 


a 


i 

6 


(5.1) 


for some e > 0, and let ( u,U ) £ A/ 1 (C2) X II be the unique solution to the forward 
problem for cr := a/|Vu|. Then 


G a {u , XJ) < G a (v, V), for all V £ II. 


(5.2) 


Moreover, if equality holds in (15.21) then (u, U) = (v, V). 

Proof. Let V £ II be arbitrary. Since (u, U ) is a global minimizer of F a as in 
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Writing 


(IA. 121) with a = a/|Vv| as shown in Theorem IA.41 we have the inequality: 

r i N r r 1 

G a (v,V) = / ci\X7v\dx + - Y' / -(v-V k ) 2 ds-2I k V k 
do. Z k=0 Ue k Z k 

= ^ J^ a W V \ dx + F v ^( V , V ) 

> \ J^ a \^v\dx + F^(u,U). 


a\S7u\dx = 


/ o 


< 


|Vu|J 

a 


| V i 


|Vt 


| Vitjcfa 


In [i IW' V “' 2<fc 


<i/ n a|V„|* + i/ o] ^|V„fe 


we also obtain 


(5.3) 


G a (u,U)= f a\S7u\dx + ^ Y^ [ —{u-U k ) 2 ds~2I k U k 
Jo 2 Ue fc 

1 /* i /» -j^ /* i 

<77 / a|Vr|dir + — / Vw| 2 dx + - / —(u - U k ) 2 ds - 2I k U k 

2 Jo 2 Jo |Vw| 2 ^ Lie fe -fe J 

= ^^«|Vw|da; + F_|_(w,t/). 


(5.4) 


From (15751) and m we conclude (15.21) . Moreover, if the equality holds in (15.21) 
then equality holds in (15.31) . and thus 


dd'a/l Vvl U) -^a/|Vf| (^5 ^0* 


(5.5) 


Since (u , [/) is a solution to the forward problem (for cr = a/|Vu|) it is also a global 
minimizer of F’ q /|vi>| over H 1 ^) x II. But (15.51) shows that (v,V) is also a global 
minimizer for F a / j vw| - Now the uniqueness of the global minimizers in Theorem IA.41 
(for a = a/\S7v\) yields (u, U) = ( v , V). 

□ 

Algorithm: We assume the magnitude a of the current density satisfies 

essinf(a) > 0. (5.6) 

Let e > 0 be the lower bound in (11.81) , and S > 0 a measure of error to be used in 
the stopping criteria. 

• Step 1: Solve (11.11 11.21 11.31) and m for cr = l, and let uq be its unique 
solution. Define 

a 1 1 

|V^’7 ’e 


g i := mm < max 
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• Step 2: For a n given: Solve POQ and (11.41) for the unique solution 

'U j n\ 

• Step 3: If 

UVUn-VUn-rllcflj) 

then define 

^ +i: =mm{max{ 1 ^- I ,e},l| (5.7) 

and repeat Step 2; 

• Else STOP. 

6. Numerical Implementations. We illustrate the theoretical results on a 
numerical simulation in two dimensions. 

6.1. An algorithm for the forward problem. Given a current pattern I £ 
II and a set of surface electrodes e k with impedances z k (taken to be constant)for 
k = 0,1,..., N satisfying m, our iterative algorithm consists in solving the forward 
problem (II. ID . (11,21) . (1 1.31) and (11.41) for an updated conductivity at each step. 

A piecewise linear (spline) approximation of the solution to the forward problem 
is sought on an uniform triangulation of the unit box [0,1] x [0,1] as shown in figure 

ICT 



1 2 3 4 


Fig. 6.1. The uniform triangulated unit box with m = 16 nodes and 18 triangles. The shaded 
triangles are adjacent to the 10 -th node. Note that the corner node 1 has one adjacent triangle, and 
nodes f., 13 and 16 have only two adjacent triangles. 


For a (square) number m of grid nodes, let T) be the set of planes supported in 
the Z-th triangle A;, for l = 1,2,..., 2 [y/rn — l) 2 . More precisely, when l = odd, 

Ti = jl - I( x - x kl )~ l(y - y kl ), ~(x - x kl ), l(y - y kl ) j , (x,y) £ A h 

where (x kl , y kl ) is the southwest grid point of the square in which A; is inscribed, and 
h is the length of the side of the square. 

When l = even, 

Ti = {l + ^-(aj-avj + jjy - y ri ), -^{y - y ri ), -^(x - x n ) \ , (x,y) £ A/, 
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where ( x ri , y ri ) is the northeast grid point of the square in which A; is inscribed, and 
h is the length of the side of the square. For example, in Figure RTT1 the triangle An 
lies in a square whose southwest grid point position is (xy, y?) and the northeast grid 
point location is (£ 12 , 2 / 12 )- 

We seek an approximation to the solution of the forward problem (ED , ot . 
OK. and (11.41) in the form 

m 

u{x,y) ~'^2u j 'ijj j {x : y), (6.1) 

3 =1 

where ipj is the sum over those planes T)’s, that are adjacent to the j-th node in the 
unit box, see figure [tTTTl By substituting m into (IA. 131) . and by selecting v = ipj , 
for j = 1 ,..., m, and V = 0, we get the set of equations 

N 

/ ctVm • Vipjdxdy + ^ — I (m — Uk)4>jds = 0, Vj = 1,..., m, (6.2) 

3 0, k—0 ~ k dek 

which is augmented with the second set of equations 

N I /• N 

-Er/ = Vj = 0,... A — 1 (6.3) 

fc=0 ° k ^ ek fe=0 


obtained by setting v = 0 and VJ? = 1 whenever k = j for k = 0,..., N — 1, and 
V ] N = -1 for j = 0,... N- 1. 

Note that forming V in this fashion is equivalent to choosing for each j = 
0,1,..., N — 1 the j-th vector for the j-th equation in (16.31) from the set 



1 


0 


0 



0 


1 


0 



... 0 


... 0 

5***5 

0 



0 


0 


1 



-1 


_ -1 _ 


-1 

> 


which is a basis for II, and then setting 

iv -1 

Un = — Uk- 

k-0 


The values of {mi, M 2 ,..., u m , Uo, U \,..., Ljv-i} are then solutions to the linear sys¬ 
tem 


A 


Ml 


0 

U 2 


0 



0 

Uo 


do — In 

Ux 


h-I N 

1 - 

1 

l_ 


1 

1 

1 

1 _ 


(6.4) 
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where the entries of A are 


N 

A i tj = / aVipj ■ Vipjdxdy + V — / tpjipids, i,j = 1,2,..., 

JQ fc—0 J e k 


m 

the entries of ’F are 
1 f 




Z N 


and 


r = 


ipids - 

Zk . 

/ i’ids, i 

' Zk 

= 1,2,.. 

., m, & k ■ 

l e o| l e ivl 

Zq ' Z N 

l e ll 


\e N -i 

Zl 


ZN-1 

l e ol 

l e i 1 I l e iV 1 

Zl ' Z N 


\eN-i 

zq 


ZN-l 

. O o 


|ejv- 2 | 

ZN- 2 

|ejv — 1 1 , 
ZN-l 


In the matrix above |e 7 1denotes the surface area of the electrode, for j = 0, ...,7V. 

For other numerical schemes for solving the forward problem we refer to E- 

6.2. Simulating the interior data. We consider a simulated planar conductiv¬ 
ity a (which models the cross section of a torso) embedded in the unit box [0,1] x [0,1]; 
see Figure [71721 on the left. The values of the conductivity range from 1.0 S/m to 
1.8 S/m. 



Magnitude of the Current Density 



Fig. 6.2. The exact conductivity with the electrode set up (left). The simulated magnitude a of 
the current density field (right). 

Two currents —Iq = I\ = 3 mA are respectively injected/extracted through the 
electrodes 

e 0 = {(x, y) e [0,1] x [0,1] : y = 0} and ei = {(x, y) £ [0,1] x [0,1] : y = 1} 

of equal impedances zq = Z\ = 8.3 mil • m 2 . 

For the given a we solve the forward problem CB, G2D, for (u,u). 

The interior data of the magnitude a of the current density field is defined by a := 
cr|V?i|; see Figure IQ1 on the right. 
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6.3. Numerical reconstruction of a simulated conductivity. Knowing the 
injected currents Iq and Ii, the electrode impedances Zq and z i, and the corresponding 
magnitude a of the current density we find an approximate minimizer of G a via the 
iterative algorithm in section [5] The iterations start with the guess uo = 1. An 
approximate solution v is computed on a 90 x 90 grid. The stopping criterion (15.71) for 
this experiment used (5 = 10 7 , and was attained with 320 iterations. An intermediate 
conductivity a v := a/|Vu| is computed using the computed minimizer v, see Figure 

m 


a : Minimal Conductivity 



0 0.2 0.4 0.6 0.8 1 


Fig. 6.3. Without any boundary measurement an intermediate conductivity can be recovered 

We note that in this example | J\ > 0.5 > 0 everywhere, thus all the level sets of 
u are connected. It follows from the arguments in Section 3 that there is a unique p 
such that u(x) = p(v(x)) and 

<r(x) = x en. 

r[v{x)) 

The function p can be determined from knowledge of u on the curve T = {(l,y) : 
0 < y < 1}, which connects the two electrodes. More precisely, for each point on T 
the function p maps the computed value of v to the measured value of u|r at the 
same point. Figure I7TT~T1 on the left shows the resulting ip on the range u(r). Figure 
16.41 on the right shows 1/p' o v, which is the scaling factor needed to obtain the true 
conductivity a. 




Fig. 6.4. Graph of the function p (left) vs. graph of 1 /(ip' o v) (right). 
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Exact Conductivity Scaled Conductivity 



0 0.5 10 0.5 1 


Fig. 6.5. Exact conductivity (left) vs. Reconstructed conductivity (right) 


In Figure [6751 the reconstructed conductivity a is shown on the right against the 
exact conductivity on the left. The L 2 error of the reconstruction is 0.04. 
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Appendix A. A minimization approach for the Complete Electrode 
Model. 

In this appendix we show solvability of the forward problem for the Complete 
Electrode Model of [25| by recasting it into a minimization problem. While this ap¬ 
proach is less general than the one given in |28| (we assume a real valued conductivity 
and positive electrode impedances), it explains how we are led to introduce the func¬ 
tional (11.131) in the solution of the inverse problem. For existence and uniqueness of 
solutions of the forward problem, the conductivity and electrode impedances need not 
be smooth: er g L°°(D) and Zk £ L°°{ek) satisfy 

essinfncr > e > 0, essinf ej , > e > 0, k = 0,..., N. (A.l) 

Let H l (£i) be the space of functions which together with their gradients lie in L 2 (D), 
and II be the hyperplane in (11.61) . We seek weak solutions to (11.11) . (11.21) (11.31) . (11.41) . 
and (11.51) in the Hilbert space H 1 ( fl) x n, endowed with the product 


((it, U), (v, V)} := f uvdx + f 
Jn J a 


V?i • 'S/vdx 


N 

E 

*;=0 


£414, 


and the induced norm 


\\(u,U)\\ := {(u, U),(u, U)) 1/2 . 


(A.2) 
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We’ll need the following variant of the Poicare inequality, suitable for the complete 
electrode model. 

Proposition A.l. Let fl c R™, n > 2, be an open, connected, bounded domain 
with Lipschitz boundary <9f 1, and II be the hyperplane in (fL?)l) . For k = 0 
let ek C dLl be disjoint subsets of the boundary of positive induced surface measure: 
|e*;| > 0. 

There exists a constant C > 0, dependent only on SI and the eu’s, such that for 
all u £ H 1 ( S7) and all U = (C/o, •••, Un) £ II, we have 

[ u 2 dx + t Ul<c( I \\7u\ 2 dx + J2 [ ( u - u k) 2 ds) . (A.3) 

k—0 V Q k—0 d ek ) 


Proof. We will show that 

In |Vu| 2 dx + Ek=o 5eS u ~ U k) 2 ds 

mf -- 

(u,u)eH H n )x u f n | X/ u \Hx + f Q u 2 dx + £f =0 


=: k > 0 . 


(A.4) 


We reason by contradiction: Assume the infimum in (IA.4I) is zero. Without loss of 
generality (else normalize to 1), there exists a sequence {{u n , U n )} in the unit sphere 
of H 1 ^) x n, || 


, £/ n )|| = 1, and such that 


0 = lim / |Vm„| 2 gU;, 

n - > °° Jn 

0= lim [ (u n — Ujf) 2 , for k = 0,..., N. 
n ^°°Je k 


(A.5) 

(A.6) 


Due to the compactness of the unit sphere in II and of the weakly compactness of the 
unit sphere in H 1 (S4) it follows that there exists some (u*,U*) £ H 1 ($7) x II with 


such that, on a subsequence (relabeled for simplicity), 

u n —>• it* iniJ 1 (D), 


U n 


U* in II, as n —> oo. 


(A.7) 

(A.8) 

(A.9) 


Since the sequence {u n } is bounded in H 1 (I4), the trace theorem implies that u ra | efc 
is (uniformly in n) bounded in H 1 ^ 2 (ek), hence also in L 1 (ek), for each k = 0, 

Using (IA.6I) and (1A.9I) in 

[ (u n - Ul) 2 ds = f 0Un - U£) 2 ds + 2(Ujf - Ul) f u n ds 
J efc J ek J ek 

+ |e fc |[(t/ fc *) 2 -(U fe ") 2 )], 


we obtain u„ 


U£ in L 2 (efe). Since 


we conclude that 


u*| efc = Uf for each k = 0, N. 
Now using (IA.5I) and (I A. 81) 


0< / |V(u n — it*)| 2 da: = / |Vu n | 2 da; — 2 / V«„ ■ Vu* + / |Vw*| 2 d:r 

Jo, J q, J q, J q, 

—> - / \\7Ui,\ 2 dx, 

Jn 


as n —> oo, 


(A.10) 
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and, since 57 is connected, 


zt* = const, in 57. (A.11) 

From (IA. 101) and (1 A. Ill) we conclude that zt* restricts to the same constant on 
each electrode, and thus Uq = U* = ... = U jy = zt*. Since U* C II, we must have 
U* = (0, ...,0) and then zt* = 0, thus contradicting (1A.7I) . 

□ 

Proposition A.2. Let 57, n, and e k C 557, k = 0, ...,7V be as in Proposition 
CO For a, and z k , k = 0,..., N satisfying (IA.1I) . and I = (Jo, ..., In) G l^ 1 , let us 
consider the quadratic functional F„ : if 1 (57) x II —» R. defined by 


N 


F a (u, U) := 


i an 


cr|Vzt| 2 da; + -^ 


k=0 1 


Zk 


N 

(zt - u k ) 2 ds - ^2 hUk■ 

k =0 


(A.12) 


Then 

(i) F a is strictly convex 

(ii) F a is Gateaux differentiable in AT 1 (57) x II, and the derivative at (zt, U) in the 
direction ( v , V”) is given by 


(DF a (u,U)\{v,V)) = f aVu ■ Vvdx + ^ / — (zt — U k ){v — Vk)ds 
JQ k—0 ^ ek 

N 

- ^ hVk (A.13) 

k =o 

(in) F a is coercive, more precisely, 

i n 

F a {u,U)> C -\\{u,U)\\--Y. I l ( A ' 14 ) 

L fc=o 

for some constant c > 0 dependent on the lower bound e in (IA.1I) . and k in (IA.4I) . 

Proof, (i) The functional has two quadratic terms, each strictly convex, and one 
linear term, hence the sum is strictly convex, (ii) The Gateaux differentiability and 
the formula (IA.13I) follow directly from the definition of F a . 

(iii) Proposition I A. 1 1 above shows that 

N 

Fffu,U)>cUu,U)\\ 2 -J2hU k , 

k—0 

where c = l|e > 0. By completing the square one obtains 

N / T \ 2 , N 

Fffu, U) > c\\u\\ 2 H1(n) +cY,(Uk-Y c ) ^ £ % 

k=0 ' / 6 k—0 

- C M 2 H1{Q) + c £ (2^“ 4^) “ 

k—0 V 7 k =0 

k—0 
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□ 

The proposition below revisits [2SJ Proposition 3.1.] and separates the role of the 
conservation of charge condition ED- This becomes important in our minimization 
approach, where we shall see that F a has a unique minimizer independently of the 
condition of ED being satisfied. However, it is only for currents satisfying ED, 
that the minimizer satisfies ED- This result does not use the reality of cr and of Zk s. 
Recall that the Gateaux derivative of DF a is given in (1A.13I) . 

Proposition A. 3. Let Q, n, ek C dfl, Zk , k = 0, and a be as in Proposi- 
tion I A. 21 

(i) If (u,U) £ H 1 (LI) x n is a weak solution to (11.11) . (11.21) . (11.31) and ED, then 
ED holds and 

(DF„(u,Uy,(v,V))=0, V (v, V) £ H 1 ^) x n. (A.15) 

(ii) If (u,U) £ H 1 ( f2) x n satisfies (1 A. 151) . then it solves (11.11) . (11.21) and ED- 
In addition, if Ik’s satisfy ED, then ED also holds. 

Proof, (i) Follows from a direct calculation and Green’s formula. 

(ii) Assume that (1A.15D holds. 

By choosing v £ Hq(Q) arbitrary and V = 0 in (IA. 151) we see that 


/ crVu ■ "Vvdx = 0. 

Jn 

Thus u £ is a weak solution of (11.11) . 

For each fixed k = 0, 1, ..., N keep V = 0 as above, but now choose v £ 
arbitrary with = 0. A straightforward calculation starting from (IA.15I) shows 

that 

/ Qu\ 

lu-U k + j y d s = 0- 

Since v\ ek were arbitrary (11.21) follows. 

Now choose V = 0 as above but v £ H 1 (H) arbitrary with v\ ek = 0 for all 
k = 0, ...,7V. It follows from (IA.15I) that 

f „pvd, = 0. 

Jon dv 

Since the trace of v is arbitrary off the electrodes ED holds. 

Finally, for an arbitrary V £ n choose v £ H 1 (Ll) with the trace v = 14 on 
each efc, k = 0,..., AT and v = 0 off the electrodes. By using the already established 
relations (11.11) . (11.21) . (11.41) and Green’s formula in (1A. 151) we obtain 





= 0 . 


On the one hand, by introducing the notation a := (ao, ...,ajv) with 


Oik ’■= 



du , 
a—ds 
ov 


k = 0,..., N, 


we just showed that a _L n. Note that so far we have not used the conservation of 
charge condition ED- 
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On the other hand, by using dm , (fT3i) . and ED in the divergence formula, we 
have 

-N n QlX f 

/ otk = o—ds = / V ■ aS/udx = 0, 

Jan ov Jn 

which yields 7$ £ II. Therefore 7$ £ II- 1 fl II = if, and (11.31) holds. 

□ 

The following result establishes existence and uniqueness of the weak solution to 
the toward CEM problem; contrast with the proof of Theorem 3.3 in [251 . 

Theorem A.4. Let LI, n, efc C dLl, Zk , for k = 0, ..., N, and a be as in Proposi¬ 
tio n A. Let F„ : H 1 (LI) x II —> R. be defined in (IA.12I) . 

(i) Then F a has a unique minimizer ( u,U) £ Ft 1 (H) x II. If, in addition, the 
injected currents Ik ’s satisfy (ED the minimizer is the weak solution of the problem 
(11.11) . (11.21) . (11.31) . and (11.41) . 

(ii) If the problem (ED , (ED. (ED,(ED has a solution, then it is a minimizer 
of F a in the whole space H 1 (LI) x II and hence unique. Moreover, the current Ik’s 
satisfy (ED- 

Proof, (i) Let 


d= inf F a (u,U ), 

and consider a minimizing sequence {( u „,, U n j} in H 1 (S2) x II, 


d < F a (u n , U n ) < d H—. 


(A.16) 


Since inf F a > — A J2k =o we h ave d ^ —oo. Following (IA.14I) . 

lim F a (u, U) = oo. 

||(u,E/)||-Kx> 


Thus the minimizing sequence must be bounded, hence weakly compact. In particular, 
for a subsequence (relabeled for simplicity) there is some (w„, U*) £ iL 1 (H) x II, such 
that 


u n — 1 u* in IL 1 (fI), and U n —> U* in II, as n oo. (A.17) 

On the other hand since F a is convex, and Gateaux differentiable at (u*,t/*) in the 
direction (u n — u*,U n — U*), we have 

F a (u n , U n ) > F a (u„ U*) + (DF a (u*, U*)-, (u n - u„ U n - U*)). (A.18) 

We take the limit as n —> oo. The weak convergence in (IA. 171) yields 
(DF a (u*, U*), (u n - u„ U n - U*)} -t 0. 

Thus d > F a (u*,U*) > d which shows that («*,!/*) is a global minimizer. Strict 
convexity of F a implies it is unique. At the minimum (u*,t/*) the Euler-Lagrange 
equations (IA. 151) are satisfied. An application of Proposition IA.3I part (ii) shows that 
(u*, U*) is a weak solution to the forward problem. 
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(ii) Proposition IA.31 part (i) shows that (u*,t/*) solves the Euler-Lagrange equa¬ 
tions, and due to the convexity it is a minimizer of F a . Due to the strict convexity of 
the functional the minimizer is unique, hence the weak solution is unique. 

□ 

Appendix B. On the regularity up to the boundary for CEM boundary 
conditions. 

If cr € C a (Cl), then interior elliptic regularity yields u £ C 1 ,a (Q). The following 
result considers the regularity up to the boundary; part a) in the proposition below was 
already proved in [3j Remark 1]. We reproduce the proof for the reader’s convenience. 
Let E := U fc=0 efc c ^ be the un i° n set °f the electrodes. 

Proposition B.l. Let f l c R" be a C 2 -domain, and let a £ C a (Ll) be C 2 - 
smooth near the boundary. Assume that the electrodes efc has Lipschitz boundary, and 
Zk £ C 2 {e k ), k = 0, ...,1V. Let ( u,U ) £ 7L 1 (S1) x II be the solution to the forward 
problem. 

Then 

a) u £ H 2 ~ e (Ll), for all e > 0. 

b) For any Xq £ dfl \ dE, there is a neighborhood Vo C Q of Xq, and a function 
vo £ H 3 ~ e (Ll) for all e > 0, such that u\v 0 = v\v 0 - 

Proof. 

a) Since u £ H 1 (Q) it follows from (11.21) that z k cr-j^u £ H 1 / 2 (E). Since 1 /{z k a) £ 

C 2 (E), we have £ H 1 // 2 (E). By [5J Theorem 11.4] the extension by zero to the 
whole boundary yields -j^u £ H 1 ^ 2 ~ e (dCl), and thus £ H 1 /, 2 ~ e (dfl). Now 

apply the elliptic regularity for the Neumann problem R5] Remark 7.2] to conclude 
u £ Ft 2 ~ e (n), for all e > 0. 

b) Let x 0 £ E, and choose r 0 > 0 be sufficiently small so that a is C 2 -smooth in 
{x £ Q : \x — xo| < 2ro} and {x £ dfl : \x — a?o| < 2r 0 } C E. 

Let Vo := {x £ S2 : \x — xo| < ro/2}. We define 

v o{x) = Xo(x)u(x), x £ fl, 

where \o £ is the cutoff function with yo(a ; ) = 1, if x £ Vo, and \o(x) = 0 if 

{x £ Ll : \x — Xo| > ro}. Then, by part a) we have 

V ■ aVvo = uV ■ aVxo + 2crVxo • Vu £ 


for all 0 < e < 1. 

Also by part a) we have that the trace u\on £ Fl 3 ^ 2 ~ e (dfl). Now, by (11.21) . 

dvo du dxo 

2 fccr— = Xo z k (T— + z k au—— 
av ov ov 

= Xo(U k -u) + z k au ^ G H 3/2 ~ e (dn), 

and thus, £ H 3 ^ 2 ~ e (dfl). Now apply the elliptic regularity [9] Theorem 7.4, 

Remark 7.2] to conclude vq £ H 3 ~ e (Ll). 

If xo £ dfl \ E, the same proof holds if we choose r o > 0 such that {x £ dLl : \x — 
Xo| < 2r 0 } D E = 0. □ 

Appendix C. The connected components of almost all level sets reach 
the boundary. Let t be one of the values for which the level set {u( x) = t} is a 
C 1 - smooth hypersurface (which is the case for a.e. t), and E be one of its connected 
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components. We show here that E D dfl ^ 0. The arguments in the proof of [22] 
Theorem 1.3] use only the interior points of f2, and thus apply to the CEM as well. 
We include them here for the convenience of the reader. 

Arguing by contradiction, assume that Ep|9fl = 0. Then E is a compact 
manifold with two connected components. Using the Alexander duality theorem in 
algebraic topology for dfMJ E (see, e.g. Theorem 27.10 in [5],) we have that R n \ 
(<9S2 [J E) is partitioned into three open connected components: (R n \ fi) [J 0\ [J 0 2 . 
Since ScSlwe have 0\ (J 0 2 = Q \ E and then dOi C dfl (J E for i = 1,2. 

We claim that at least one of the dO\ or dOi is in E. Assume not, i.e. for 
each i = 1,2, dOif] |<9fl ^ 0. Since dfl is connected (by assumption), we have that 
O 1 UO 2 U is connected which implies 0 \ (J O 2 (J ( R n \ U) is also connected. By 
applying once again Alexander’s duality theorem for E C R n , we have that R n \ 
E has exactly two open connected components, one of which is unbounded: R n \ 
E = OooU^o- Since OiU^zLK^" \ ' s connected and unbounded, we have 
0\ (J O 2 (J{R n \ U) c Ooo , which leaves Oq C R n \ (Oi (J O 2 U(^ n \ ^)) C E. This 
is impossible since Oq is open and E is a hypersurface. Therefore either 0\ or O 2 or 
both has the boundary in E. 

To fix ideas, consider dO\ C E. If this were the case, then we claim that u = t 
in 0\. Indeed, since 0\ is an extension domain [dO\ has a unit normal everywhere) 
the new map u defined by 



u(x), x € Q\0 1 , 
t, x £ 0 1 , 


is in H 1 (U) f) C(fl) and decreases the functional (11.1311 . thus contradicting the mini¬ 
mizing property of u. Therefore u = t in Oi, which makes |Vit| = 0 in 0\. Again we 
reach a contradiction since the set of critical points of u is negligible. 


These contradictions followed from the assumption that E f) dfl = 0, and there¬ 
fore E f) dfl ^ 0 
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